Optimization of the current pulse for spin-torque switches 
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We address optimization of the spin current intensity profile needed to achieve spin torque switch- 
ing of a nanomagnet. For systems with Ohmic dissipation we prove that the optimal current drives 
the magnetization along the trajectory, which is exact time-reversed replica of the relaxation trajec- 
tory towards the equilibrium. In practice it means that the optimal current is very nearly twice the 
minimal critical current needed to switch the magnet. Pulse duration of such an optimal current is 
a slow logarithmic function of temperature and the required probability of switching. 

PACS numbers: 75.70.-i, 85.75.-d, 75.75.Jn 



The spin torque effect, proposed by Slonczewski[l| and 
Bergerpj, is a subject of intense study 0-[l2|> because 
of its ability to cause magnetic switching in ferromag- 
netic structures. This opens the possibility of using it 
for high capacity and low volatility data storage appli- 
cations. One of the key problems is minimization of the 
energy dissipation during the switching event. Most of 
the research centers around optimizing material and ge- 
ometric parameters of the nanomagnetic structures how- 
ever some look at how the shape of the spin current pulse 
can be used 
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Here we focus on a much simpler 
aspect: optimizing intensity and time dependence of the 
spin current pulse. A very high current density, while 
achieving fast switching, results in large Joule heat. On 
the other hand, the low current density requires a long 
switching time, which again brings substantial heating. 
Clearly there is an optimal switching protocol, with min- 
imal deposited heat. 

Surprisingly there is a generic prescription for such an 
optimal current protocol, which optimizes Joule losses ir- 
respectively to specific parameters of the structure. We 
show here that the optimal spin current is such that the 
sample magnetization retraces its relaxation trajectory 
in reversed time. The relaxation trajectory is the magne- 
tization history starting from the unstable energy maxi- 
mum and relaxing to the stable equilibrium in the absence 
of any external spin current. In practice it implies that 
the optimal current intensity is almost time-independent, 
which is about twice the minimal critical switching cur- 
rent. The duration of such an optimal current pulse is 
non-universal and is sensitive to the confidence level of 
the switching probability as well as temperature. 

To be specific we focus on a mono-domain soft ferro- 
magnetic layer, whose magnetization, M(i), dynamics is 
described by the Landau-Lifshitz-Gilbert equation 
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where Mcons corresponds to conservative motion along 
the Stoner-Wohlfarth (SW) orbitals. The dissipative 
Gilbert term Mdiss is perpendicular to SW orbitals 
and causes the magnetization to relax towards the easy 



axis/external field equilibrium direction. The last term, 
Mst , describes the spin torque effect and has components 
both perpendicular and parallel to SW orbitals. These 
three torques may be written as 

Mcon. = X He//, 

Md^.s = -7am X (M X He//) , (2) 

Mst = ~^J{t)M X (m X z) . 

Here He// = —VmE/iiq is the effective magnetic field 
of the system and E is the energy given by 
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Mo 2M, 2Ms 
where M — mMj and Ms is the saturation magnetiza- 
tion, is the easy axis anisotropy field strength with 
easy axis along the z direction and iJ^ is the strength of 
the easy z — y plane anisotropy field. The Gilbert damp- 
ing constant is a and J{t) represents the strength of the 
spin current, which is polarized along the z direction. 
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FIG. 1. Azimuthal angle 6 vs. time (full line) and energy vs. 
time (dashed line) for magnetization reversal of the system 
with Q = 0.01, HI = 0.028M,, HI = Ms, and J = O.OIM,. 
Time is measured in units of to = (^Ms)^^ , while the energy 
in units of Eo/n — ^oHlMs/2n. 

To analyze the system's dynamics it is convenient to 
exploit the time scales separation between the preces- 
sional time and the switching time. Figure [1] shows the 
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azimuthal angle evolution during the switching process. 
The system goes through many revolutions as it moves 
out of its initial 9 = direction towards the switching 
point at 9 = tt/2, and further during its relaxation to- 
wards 9 = IT. At the same time the energy E{t), Eq. ([3]), 
is a rather smooth function of time, which first increases 
(thanks to spin-torque) and then relaxes to the new min- 
imum. This observation suggests to integrate out the 
fast degree of freedom by averaging over each SW or- 



bital, labeled by its energy E, Ref. [15[. This procedure 
results in the effective Id dynamics, with energy as the 
only "coordinate" 
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lM,P{E) 
-aU,ss{E) + J{t)rh-lst{E), 



(4) 



where P{E) is the period of the SW orbital with energy 
E, while Idiss{E) and lst{E) are the contour integrals 
along this orbital of the dissipative and spin torques, cor- 
respondingly 



lMss{E) = ^^j> [dMxHe//] 
lst{E) 



MsP{E) 



dM X M. 
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FIG. 2. (Color online) Energy flow E, Eq. (|4]), in arbitrary 
units as a function of energy E for various strengths of the 
spin current J for the system with a = 0.01, = 0.028Afs 
and = Ms. From top to bottom J = 0.012Af^ (red), 
Jc = 0.0066Ms (orange), J = 0.0052M., (green), J = 
(blue). Black dashed line indicates the optimal spin current 
J opt, Eq. ((S)), as a function of energy in arbitrary units. 

Figured] shows how E behaves under various strengths 
of the applied spin current J. Without the spin cur- 
rent the energy flow is of purely dissipative nature 
—aIdiss{E) < 0, which forces the magnetization to relax 
towards = 0. On the other hand, for the spin current 
above a certain critical value J > Jc, the energy flow is 
positive, resulting in an eventual magnetization switch. 



as shown in Fig. [TJ There is also a narrow interval of the 
spin currents 0.82Jc ^ J ^ Jc, where the energy flow is 
initially positive and changes sign at some energy below 
the maximum. This leads to stable precession around a 
proper SW orbital observed in a number of studies 

In order to find the optimal switching protocol J{t) we 
require that the Joule heat / dtJ^{t)R, deposited dur- 
ing the switching event, is minimal. Here we assume for 
simplicity that the total current is proportional to the 
spin-current J. We also assume that the resistance R of 
the device stays approximately constant as the magneti- 
zation direction evolves. To ensure the minimum of the 
Joule heat one needs to require that the ratio of the en- 
ergy gain E{t) to the energy loss J^{t)R is maximized at 
each instance of time. Substituting the energy equation 
of motion @ and taking derivative with respect to the 
current, one finds the following condition 
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Solving for the spin current J, yields 

Idiss {E) 



Jopt = 2a 



m-ht{E) 



(7) 



(8) 



Therefore the optimal current is exactly twice the one 
needed to nullify the energy flow (U). Substituting this 
back into equation Q , one finds that the latter takes the 
form 



E : 



-aIdiss{E) 



(9) 



This is exactly the same as in Eq. Q without the ex- 
ternal spin current J, but with the time being reversed. 
Therefore the optimal current protocol is such that it re- 
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FIG. 3. (Color online) Joule heating during the switching pro- 
cess in arbitrary units vs. spin current for various confidence 
levels p. From top to bottom p = 0.99 (red), p = 0.5 (green), 
p = 0.2 (blue). Here a = 0.01, = 0.028Ms, Hj^ = Ms and 
Eo = UOkbT. 
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verses the purely relaxational trajectory in energy of an 
isolated system. Such a relaxation trajectory must be 
thought of as starting at the energy maximum (i.e. the 
switching point) and winding down towards the minimum 
at 9 = 0. Since the ratio of the dissipative and the spin- 
torque currents along such a trajectory is roughly energy 
independent, see Fig.[2l equation ([8]) implies that the op- 
timal current is very nearly a constant given by twice the 
critical current Jc- This is shown in Fig. [3] via Monte- 
Carlo simulations of the Joule heating as a function of 
spin current for various confidence levels. Thermal fluc- 
tuations enter as a random component Urand of the ex- 
ternal magnetic field with zero average and mean square 
value (H.land) — '^jMsakBTAt/ Ho [20| (we assume tem- 



scales as 



peratures large enough to ignore shot noise [ll|). At low 



confidence levels J opt however at confidence levels 

nearing 100% we see the predicted Jopt ~ 2Jc. 




FIG. 4. Simulated optimal switching time as a function of the 
confidence level (crosses), Eq. (|10p (full line). Same parame- 
ters as in Fig. [3] with J = Jopt- 



While we have established that the optimal current 
pulse is Jopt ~ 2Jc, we have not discussed yet its dura- 
tion. Here we run into a problem. Indeed, the dissipative 
relaxation trajectory formally takes an infinite time to 
approach the energy minimum a,i — Q. Thus its time- 
reversed, the optimal trajectory, is infinitely long too. 
What saves the day is the fact that the initial direction 
of the magnetization is not exactly aligned along = 0. 
Instead it is scattered around this direction according 
to the equilibrium thermal distribution cx e"^*"'/'^^-^, 
where Eini is the starting energy counted from its min- 
imal value at 6' = 0. Thus the times needed to reach 
the switching point, which are oc ln{Eo / Eini) where 
Eq = fioH^Ms/2 is the hight of the energy barrier, are 
also scattered accordingly. As a result, any pulse of a 
finite duration achieves switching only with a certain 
probability p < 1. If this probability, i.e. error toler- 
ance, is specified such that 1 — p ^ 1, all realizations 
with Eini > ksT^l — p) should undergo the switch. This 
dictates that the duration of the optimal current pulse 
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a \kBT{l-p) 



(10) 



where tg = {jMs)~ . This result is shown in Fig. |4]for 
higher switching probabilities and is in a good agreement 
with the simulated data. 

The statement that the time-reversed trajectory opti- 
mizes the switching process is not restricted to the aver- 
aging procedure employed in this paper. One can imag- 
ine sufficiently complex switching protocols, utilizing the 
precessional freque ncy to achieve such an optimal path. 
For example, Ref. |13j used current pulses resonant with 
the orbital period to excite the system. In the realm 
of relatively low frequency protocols however, twice the 
critical current pulse provides the best approximation to 
the optimal strategy. The duration of this pulse weakly 
depends on the confidence level and can be estimated 
according to Eq. (fTO| . 
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